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Abstract. The article describes an investigation of the effectiveness of genetic al- 
gorithms for multi-objective combinatorial optimization (MOCO) by presenting an 
application for the vehicle routing problem with soft time windows. The work is 
motivated by the question, if and how the problem structure influences the effec- 
tiveness of different configurations of the genetic algorithm. Computational results 
are presented for different classes of vehicle routing problems, varying in their cov- 
erage with time windows, time window size, distribution and number of customers. 
The results are compared with a simple, but effective local search approach for 
multi-objective combinatorial optimization problems. 



1 Introduction 

The requirements of modern logistics are widespread. It seems natural that 
for practical problems several objectives have to be optimized. The quality 
of service and products and the corresponding cost are only a small subset 
of necessary objectives that should be considered [25] , 

The distribution of products, a part of the supply chain in logistics, can be 
modelled as a vehicle routing problem (VRP). The objective is to find a cost- 
minimizing set of routes from a depot to serve a number of customers with 
known demands. Each customer is serviced by exactly one of the vehicles 
which capacities are not allowed to be exceeded. 

The model of the vehicle routing problem with time windows (VRPTW) gen- 
eralizes the basic problem by introducing time windows for each customer, 
defining an earliest and a latest possible time of service. In the case of the ve- 
hicle routing problem with soft time windows (VRPSTW), violation of these 
time windows is possible and does not lead to infcasibility of the solution 
PQ. In our approach, the possible violations of the time windows are mod- 
elled by introducing objective functions, leading to a multi-objective model 
that is able to describe the requirements of practical problems more detailed 
compared to single-criterion optimization approaches. The goal is to find all 
solutions minimizing the defined objective functions. As these functions are 
often conflicting, the concept of Pareto optimality is used to determine ap- 
propriate solutions. 

The VRPTW [3T], as well as the VRP [23], has been shown to be NP-hard in 
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terms of complexity. In combination with the practical relevance, this is one 
of the reasons why there is still a lot of ongoing research. Earlier approaches 
include heuristics and optimization algorithms. Heuristics are focused on tour 
construction [31] , tour improvement [25] or both in parallel |26l27j . 
Exact optimization methods successfully apply branch-and-bound methods 
[22] or integer linear programming based on Lagrangian relaxations [11|12|13] . 
However, only smaller problem instances up to 50 customers can be solved 
reliable [21]. 

More recent, metaheuristic strategies have been applied to the VRP and 
the VRPTW. They include Simulated Annealing [5124] . Tabu Search [30134] . 
Genetic Algorithms (GAs) [35136] . Ant Colony Optimization [4], hybrid ap- 
proaches [37 38] and newer local search concepts like e.g. Guided Local Search 

M- 

Nevertheless, the focus of the treated models is in the observed cases on 
minimizing a single objective, the cost of the solution, by minimizing the dis- 
tances travelled by the used vehicles, they do not address a more practical, 
multi-objective formulation with a relaxation of the restrictive time windows. 
On the other side, there is an increasing interest on applying evolutionary 
based optimization techniques to multi-objective problems. Evolutionary al- 
gorithms are regarded to be well suited for multi-objective problems since a 
set of alternatives is used in the optimization process and the goal is the ap- 
proximation of a set of Pareto optimal solutions. The idea is the convergence 
of the whole population towards the efficient frontier. 

Since a first sketch of an idea by Goldberg [17] , most of the present work on 
multi-objective optimization using evolutionary algorithms are focused on the 
integration of the objective functions in the calculation of the fitness values of 
the solutions, ranging from scalarizing functions to Pareto-based techniques. 
Here to mention are VEGA of Schaffer [29], MOGA of Fonseca and Fleming 
PU, NPGA of Horn et al. [191, NSGA of Srinivas and Deb [32] and SPEA 
of Zitzler [41] . These approaches include the ideas of fitness functions/fitness 
sharing, niching, mating restrictions and elitism to define algorithms that 
maintain diversity within the population, overcome local optima and finally 
converge towards the set of efficient alternatives. For a detailed overview 
concerning these topics, the interested reader is referred to Coello and Van 
Veldhuizen |7|39] . Although big progresses were made and a huge variety of 
applications were already successfully developed [5], an application to the 
vehicle routing problem under multiple objectives is still missing. To define 
an appropriate algorithm, the interchange between the problem structure, 
the configuration and the behavior of the algorithm has to be studied. Espe- 
cially for the correct use of crossover operators in multi-objective optimization 
problems, computational results are only for a few examples available [3]. 
This paper fills the gap in the described field of research. A practical, multi- 
objective model of a VRPSTW is presented and a genetic algorithm is defined 
to solve different problem instances, varying in their time window coverage, 
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time window size, number and distribution of customers. Several crossover 
operators are tested and computational results are given for the different in- 
stances. To compare the results obtained by the GA, a simple local search 
approach for multi-objective combinatorial optimization is presented. 



2 The multi-objective vehicle routing problem with 
soft time windows 

2.1 Multi-objective optimization 

The goal of a multi-objective optimization is to 



R G Q is a solution of the problem and belongs to the set of all feasible 
solutions fi. The objective functions gi{R) 1 ...gk{R) map the decision R into 
the objective space, leading to a objective vector G(R). As often conflicting 
objective functions gk{R) are considered, minimization does not lead to a 
single optimal solution but is understood in the sense of efficiency (or Pareto 
optimality) [39] . 

Definition 1. An objective vector G{R) is said to dominate G(R'), \igi(R) < 
gi(R')Vi e {1, k} A3i e {1, k}\gi(R) < gi(R'). We denote the domina- 
tion of a vector G(R) to the vector G(R') with G(R) -< G(R'). 

Definition 2. A solution R G fl is said to be efficient or Pareto optimal, if 
SR 1 G Q\R' -< R. Set set of all solutions fulfilling this property is called the 
Pareto set P. 

From the description of the multi-objective optimization problem in equa- 
tion [1] we derive in combination with the definitions [1] and [3] the goal to find 
all R G P. 

2.2 Problem description 

The vehicle routing problem with soft time windows can be described as fol- 
lows: 

A known number of customers have to be delivered from a depot with a 
known amount of goods for which an unlimited number of homogeneous ve- 
hicles is available. It is assumed that each customer is visited by exactly one 
vehicle and a loading and a travelling constraint exists for the vehicles. A 
soft time window is associated with each customer, defining a desired earliest 
and a latest time of service. Violation of these time windows does not lead 
to infcasibility of the solution. With respect to the soft nature of the time 
windows, it is assumed that service is done immediately after the arrival of 




(1) 
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the vehicle. 

The objective of the problem is to maximize quality of service and to minimize 
cost, such that the requirements of the customers and the side-constraints are 
met. It is obvious, that the violation of the time windows has to be minimized 
in order to achieve a high quality of service. This can be done by minimizing 
the number of time window violations and the time window violations itself, 
measured in time units. 

The cost consist of a fixed part, induced by the number of used vehicles and 
a variable part, caused by the route length and the travel time. 



2.3 Notation of the model 

The customers and the depot are modelled as vertices V = {v , v\, wjv} in 
a graph G ~ (V, A), connected by a set of arcs A. In V, the depot is indexed 
by and the total number of customers is N. 

With each arc (vi , Vj ) from node Vi to Vj , a nonncgative travel time t Vi iVj , 

a distance d VuVj and a cost c VuVj is associated. It is assumed that t Vi>Vj = 

d Vi , Vj = c ViiVj V(vi,Vj) G A. 

u Vi : unloading time at customer Vi,Vi E V. 

d Vi : demand of customer t>i,Vi E V. 

d max : capacity constraint of the vehicles. 

]: time window at customer Vi,Vi G V. a Vi defines the desired earliest, 
b Vi the latest time of service. It is necessary that b Vi > a Vi Vvi G V. 
t%.: arrival time of the vehicle at customer Vi,Vi E V. 

Tk = [[l]fc, [2]fc, [n/c]/c]: route r^, representing customers in a given se- 
quence. To assure, that each customer is included in exactly one route, it is 
necessary that n ry — Vfc ^ k' and | Urj. | = | V \ {v } |. 
R = {n, r m }: set of routes of the solution. The number of routes m can 
be greater than one and does not have to be equal in every solution. 
The time t(rk) to travel route can now be obtained by calculating 



As the time window [a Vo ,b Vo ] defines the interval in which the vehicles are 
available, it is necessary that 




t(r k ) < b Vo Vr fe G R 



(3) 



The quantity q(rk) of transported goods on route is given by 




(4) 



i=i 
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The capacity constraint must be met, so 

q{r k )<d max Vr fe Gi? (5) 
w(Vi) measures the time window violation of customer Vi,Vi G V. 

w(vi) = max(max(0;a„, — i". ); max(0; t%. — b Vi )j (6) 

The fact, that the time window of customer Vi is violated is indicated by 

u(vi): 

t .\ f 1 : w(i) > 
«(») = n . .1\ - n ( 7 ) 



: w(i) = 
The objective functions are given by 



\R\ 

gi (R) = J2t(r k ) (8) 



k=l 



g 2 (R) =| R | (9) 



N 



g 3 (R) = ^2w(Vi) (10) 



2=1 

N 



g 4 (R) =J2<Vi) (11) 

i=l 

The objective functions © and ((9]) minimize the costs associated with the 
solution. Function (jHJ) minimizes the length of the routes, travelled by the 
vehicles. In addition to that, function @ minimizes the number of used 
vehicles. To maximize the provided service it is desired to minimize occurring 
time window violations. The objective functions (fT0|) and (fTTj) express this 
circumstance by the measured deviations in time units and by the number 
of violated time windows. 

According to the definitions (JTJ) and ^ given in section 12.11 the goal is to 
find all Re P. 
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3 A genetic approach 

3.1 Encoding technique, initialization 

We use a string of N genes to encode the solutions of the problem. Each gene 
represents a customer of the problem and the sequence of them corresponds 
to the sequence of visiting the customers. A possible string, e.g. "52137 
9 8 4 6 10" is then partitioned to a set of tours, e.g. n = [v$, «2, v%], r% = 
[i>3, v-j, vg, vg], ?"3 = [«4, Vq, Vio], assuring feasibility with respect to the defined 
side constraints ([3]) and (O in section [2~3l The decoding technique inserts the 
next stored customer of the representation in the current tour as long as the 
side constraints are not being violated. 

The initial population inherits n pop individuals i, consisting of randomly gen- 
erated strings. To assure a sufficient potential of genetic information and avoid 
premature convergence, n pop is chosen equal to 500. 



3.2 Fitness function 

Most investigations on genetic algorithms for multi-objective optimization 
are based on examining the effectiveness of certain approaches of fitness func- 
tions. As we have a different focus, we restrict ourselves to the use of a single 
function, similar to the the well known fitness function proposed by Fonseca 
and Fleming [Ti] , 

For individual i of the current population, the number of individuals whose 
corresponding objective vectors dominate the objective vector of i is indi- 
cated by £j. All currently not being dominated alternatives receive & = 0. 
The fitness f(i) of individual i is derived using a linear transformation with 
£max _ max ^^ an d two external given parameters f mm and f max ; 

i 

f(i) = f max - 5i!L_ i 1 (12) 

The parameters f min and f max are set to f min = 1 and f rnax = 100. 



3.3 Genetic recombination operators 

As we want to study the effectiveness of different crossover operators, several 
of them are implemented and tested. From the first approaches of genetic 
crossover operators, improvements were made by incorporating problem spe- 
cific knowledge into the recombination process. This is in the case of the 
travelling salesman problem as well as in the case of the vehicle routing prob- 
lem e.g. done by considering the distances of the customers to each other and 
selecting between two possible genes the route minimizing part. 
As we deal with more than one objective, problem specific orientation during 
the crossover process is not as obvious. We therefore restrict the investigation 
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in this first approach on problem independent working operators. 
Among the best studied crossover operators for permutation based encodings 
are the "partially mapped crossover" (PMX) [16] , the "order based crossover" 
(OBX) [33J and the "uniform order based crossover" (UOBX) [5]. We apply 
them with a crossover probability of p cross = 1, 

Crossover operators introduce implicit mutations to the offspring. To distin- 
guish between the true results obtained by crossover, no mutation operator 
is used in the basic test runs. However, for comparison reason a simple "swap 
mutation" exchanging two randomly picked genes is examined in a separate 
test run in combination with the UOBX as this crossover operator is being 
described as the best among the problem independent operators for the con- 
sidered problem type. In this configuration (UOBXA2EX) , it is applied with 
a probability ofp mut = 0.1 for the whole individual. Thus, the probability of 
an exchange of a single gene is 2p N . 



3.4 Elitism 

It has been shown, that elitism plays in general an important role in genetic 
algorithms and in particular in multi-objective GAs, see e.g. [41] . Various 
approaches exist to implement an elitist strategy, ranging from preserving 
the k best individuals in the population to storing the best individuals in an 
external set. 

We use a simple, but effective strategy to obtain an elitist genetic algorithm 
by using overlapping populations. This principle has been introduced as the 
so called "steady state" genetic algorithm. From each generation to the next, 
the population is preserved and new individuals are only inserted if they 
improve the average quality of the population. Duplicates are not allowed, 
including duplicates of encoded strings as well as duplicates of the corre- 
sponding solutions in the decision space. The improvement of the population 
is in our case measured by the ^-values. This means that the complexity of 
the necessary calculations to obtain the ^-values by pairwise comparison is 
2 n P°P. 



3.5 Termination criterion 

The algorithm can terminate after a defined time, a given number of calcu- 
lations or after the algorithm has converged to a given goal and does not 
produce any improvements for a defined time or number of calculations. It is 
certainly interesting to perform calculations until a convergence to a local or 
global optimal front is reached. The termination criterion is chosen according 
to this goal. The genetic algorithm stops after 10,000 iterations with no found 
new individual i having a £ j = 0. 
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4 A local search approach 



A simple local search approach is used to compare the results obtained by 
the genetic algorithm. It is based on the encoding scheme of section [3Tl and 
can be described as follows: 



1 initialize 

1.1 generate starting solution R 

1.2 encode starting solution R 

1.3 set pwvrox = { R y 

2 repeat until neighborhood of R € p a PP rox ± s investigated 

\/J£ g papprox 

2.1 select new solution R 6 p a PP rox 

2.2 generate neighborhood of R 

2.3 update p a PP r ° x 

2.4 set neighborhood of R as investigated 



The algorithm starts from a random generated solution. The calculation of 
the initial solution is done as in the initialization procedure of the genetic 
algorithm to assure proper comparison results. A set p a PP rox is used to store 
the set of currently not dominated solutions. According to the given defini- 
tions ((TJ and (|2|) in section 12.11 it has to be noticed that these solutions do 
not necessarily fulfill the requirements of Pareto optimality. However, as we 
do not know the real efficient solutions and the true Pareto set p true t we 
consider the nondominated solutions as an approximation of the Pareto set 

ptrue 

The algorithm then starts in step 2 to investigate the neighborhood of the 
solutions in p a PP rox . A solution R is randomly selected from the set p a PP rox 
and the neighboring solutions are calculated. Procedure 2 . 3 then inserts new 
nondominated solutions of the neighborhood of R into the set p a PP rox and 
removes dominated ones. 

The local search procedure terminates when the neighborhood all solutions 
R G p a PP rox has been investigated and as result a set of locally optimal so- 
lutions with respect to the used neighborhood is obtained. 
An approach based on the neighborhood proposed by Lin is used to per- 
form local search. Every possible substring of the encoding is extracted, re- 
versed and reinserted, resulting in N ( N ~ 1 ') calculations. The basic ideas of 
this method come from the local search descent method, incorporating multi- 
objective storage of nondominated solutions. We call this approach therefore 
multiple-objective local search descent (MOLSD). 
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5 Computational experiments 

5.1 Test problems 

To investigate the influence of the problem structure on the effectiveness 
of different configurations of the algorithms proposed above, a set of test 
instances incorporating different characteristics of the problem is necessary. 
For that purpose, 40 test instances based on the data given from Solomon [3T] 
were created. Each instance can be described using an a; (3; 7; ^-classification 
scheme, representing the following configurations: 

• a: Distribution of the customers. 

We consider test instances with clustered data sets (a = C) and with 
random distributed customers (a = R). 

• (3: Number of customers. 

Problems with 20 and 30 customers were defined. 

• 7: Coverage with time windows. 

This attribute measures the relation of the number of customers having 
time windows to the number of customer without time windows. 

• 5: Average time window size. 

The size of the time window at customer Vi is defined as by- Q"Vi • 

Example: The string "C;20;0.70;60" describes the instance with a clustered 
distribution of 20 customers, having a density of 70% with time windows and 
an average time window size of 60 time units. 

We expect the problems with random distribution of the customers and a high 
density of small time windows to be more complex to solve than problems 
with less and bigger time windows. This should affect the speed of convergence 
of the algorithms and the quality of the obtained solutions. 
The neighborhood of the defined problems with 30 customers contains 2.29 
times more alternatives than the one of the problems with 20 customers. An 
increase of the necessary evaluations of at least this factor is being expected. 

5.2 Measures of effectiveness 

Each of the performed test runs leads to an approximation p a PP rox of the set 
of efficient alternatives P. Due to the complexity of the problem, we do not 
know the true Pareto set P. However, to obtain a measure of effectiveness, we 
compute a reference set P re f including all the found nondominated solutions 
during the test runs. The effectiveness of the approximations is then oriented 
on the approximation of this reference set P re f . 

To determine the distance between two sets of alternatives, a distance func- 
tion measuring the distance of the underlying elements according to the the 
proposal of Czyzak and Jaszkiewicz [5] is used. The distance measure can be 
described as a weighted achievement scalarizing function, using it)j = -g- as a 
weight for the goals based on the spread Aj of the objective function values 
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of objective j in the reference set P re ' . For an alternative x £ p a PP rox to an 
alternative y G P re f , the distance is measured by equation (fT3| . 

c(x,y) =^max fe |0;u; i (^(a;) (13) 

Two approaches are used to obtain the distance of the approximation p a PP rox 
to the reference set P re f . We measure the average distance by equation (|14p 
and the worst-case distance by equation (|15jl . 

di = j-r min \c(x,y)\ (14) 

pre/ xeM I ^ ,y 'S V } 

c?2 = max { min{c(x, y} } (15) 

As our goal is to minimize the distance between the obtained approximations 
papprox anc j fag reference set P re ' , the distances d\ and di are to be mini- 
mized. 

To come to stable and reliable conclusions, average values of d\ and d% of 
several test runs with the same configuration are computed. 



5.3 Results 

The test runs were performed on a personal computer, equipped with an 1 
GHz Intel Pentium III Processor and 384 MB RAM. Each algorithmic con- 
figuration was executed 100 times for each test problem in order to obtain 
stable results. The average number of evaluations of the test runs and the 
performance of the implemented heuristics is given in table[T]in the appendix. 
Resulting from the lower computational effort, the MOLSD is able to evalu- 
ate about 4.70 times more solutions per second than the GAs. Each second, 
1,385 evaluations were computed by the MOLSD. 

As expected, all algorithms tend to converge slower in instances with dense 
and small time windows and random distributed customers. It is worth to 
mention, that the increase of the necessary evaluations of the MOLSD from 
a (3 — 20 to a (3 = 30 was about factor 9.14, a huge increase compared to 
the additional computations of the GAs, having factors of 1.80 for the PMX, 
1.96 for the OBX, 1.78 for the UOBX and 2.63 for the UOBXA2EX. 
The results, given in table [2] and [3] in the appendix indicate a strong depen- 
dency of the approximation quality on the distribution of the customers a. 
The examined crossover operators PMX, OBX, UOBX itself as well as the 
genetic algorithm including the swap mutation UOBXA2EX were in most 
instances with an a — C not able to outperform the simple local search 
approach MOLSD. Both defined measures d\ and di are affected by this 
property and the distance between the values of the MOLSD to the values of 
the GAs is increasing with falling 7 or rising S. We conclude, that VRPSTW 
with a structured distribution of the customers are more easy to be solved 



Genetic Crossover Operators 201 



with simple local search while the examined configurations of the genetic 
algorithm tend to get stucked in local optima. This effect grows with a de- 
clining complexity of the underlying problem, expressed by a low coverage 
with time windows 7 or wide time windows 6. Similar results have already 
been reported for the single objective VRPTW [2] and can be confirmed for 
the multi-objective formulation. 

For VRPSTW with a random distribution of the customer we observe the 
contrary. In most test instances, the best results for the average distance d\ 
are obtained by the UOBXA2EX. However, the advantage of the UOBXA2EX 
to the MOLSD is decreasing with falling 7 or increasing 5. For some problems 
with large time windows S and a small coverage 7, e.g. the "R;30;1.00;80/95/115" , 
MOLSD produces better results. 

The results of d,2 indicate, that the average worst-case deviation of the ap- 
proximation p a PP rox to the reference set P re f is superior in the approach of 
the MOLSD, except for some smaller instances with (3 — 20 and the hard 
"R;30;1.00;10". 

A closer investigation of the results shows, that the superior results of the 
UOBXA2EX is due to the use of the swap mutation operator. Comparing the 
results for the crossover-only GAs, it can be seen that with respect to the 
distances d\ and d-i OBX and UOBX perform in most instances better than 
PMX. On the other hand, a clear distinction between OBX and UOBX can't 
be concluded. 

The results show a strong influence of the mutation operator on the effec- 
tiveness for the genetic algorithm. Although the GAs lead to weaker results 
in clustered data sets than the proposed local search approach, this behavior 
seems to be typical for the presented multi-objective optimization problem. 

6 Conclusions 

The goal of the paper was to investigate the influence of the problem struc- 
ture on the effectiveness of different configurations of a genetic algorithm. 
A multi-objective formulation of the vehicle routing problem with soft time 
windows has been presented, describing the requirements of practical logis- 
tics more detailed than single criterion models. For comparison reasons, an 
approach of a local search approach incorporating storage of multiple non- 
dominated solutions was proposed. 

The results show, that the problem structure, especially the distribution of 
the customers a, is crucial for the behavior of the investigated algorithms. 
The genetic algorithms outperform the local search approach especially well 
in complex problems with dense and small time windows and a random dis- 
tribution of the customers. However, the influence of the mutation operator 
seems to be stronger than expected and described in other publications on 
GAs. 

An important role of simple (local search based) mutations can indicate, that 
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the studied crossover operators themselves are comparable weak for the multi- 
objective formulation of the problem as they do not recombine the desirable 
structures of the underlying model. Nevertheless, specific formulations of par- 
ticular multi-objective operators are still missing. A combination of genetic 
operators with local search heuristics is consequently a logical conclusion of 
the obtained results. 
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A Appendix 



Table 1. Average number of evaluations. 

a;f3;j;8 MOLSD PMX OBX UOBX UOBXA2EX 

C;20;1.00;60 38,730 108,750 54,375 75,625 80,000 

C;20;0.70;60 43,643 98,125 80,625 98,125 95,625 

C;20;0.45;60 35,488 103,750 82,500 89,375 96,250 

C;20;0.30;60 26,938 71,875 56,875 64,375 76,875 

C;20;1.00;120 36,814 105,625 50,625 66,875 48,125 

C;20;1.00;180 42,047 175,625 63,750 45,000 138,750 

C;20;1.00;240 26,125 73,750 55,625 46,875 58,125 

C;20;1.00;360 23,978 74,375 36,250 35,000 34,375 

R;20;1.00;10 90,060 425,625 381,875 438,125 474,375 

R;20;0.70;10 87,199 385,625 398,125 353,125 411,875 

R;20;0.45;10 90,011 386,875 301,875 303,750 433,125 

R;20;0.30;10 46,364 196,250 204,375 238,750 231,875 

R;20;1.00;30 90,011 440,625 350,625 404,375 426,875 

R;20;0.70;30 74,997 361,250 251,875 345,000 367,500 

R;20;0.45;30 59,113 211,875 201,875 178,125 266,250 

R;20;0.30;30 22,067 106,250 98,125 115,000 142,500 

R;20;1.00;60 61,191 315,625 260,000 259,375 278,750 

R;20;1.00;80 44,650 233,125 171,875 225,000 226,250 

R;20;1.00;95 51,897 231,875 206,250 196,875 286,875 

R;20;1.00;115 24,271 144,375 118,125 115,625 145,000 

C;30;1.00;60 510,995 177,500 155,000 185,000 301,250 

C;30;0.70;60 367,880 156,250 141,250 127,500 220,000 

C;30;0.45;60 372,534 220,000 261,250 240,000 541,250 

C;30;0.30;60 211,845 147,500 178,750 160,000 230,000 

C;30;1.00;120 515,736 143,750 156,250 203,750 180,000 

C;30;1.00;180 532,005 271,250 120,000 146,250 246,250 

C;30;1.00;240 335,864 188,750 82,500 82,500 186,250 

C;30;1.00;360 253,083 128,750 71,250 73,750 101,250 

R;30;1.00;10 754,551 958,750 466,250 585,000 955,000 

R;30;0.70;10 642,278 650,000 503,750 453,750 862,500 

R;30;0.45;10 550,319 673,750 392,500 460,000 1,383,750 

R;30;0.30;10 276,878 391,250 435,000 315,000 403,750 

R;30;1.00;30 809,535 666,250 445,000 321,250 1,025,000 

R;30;0.70;30 583,814 658,750 478,750 326,250 1,023,750 

R;30;0.45;30 471,105 405,000 381,250 256,250 503,750 

R;30;0.30;30 155,643 236,250 246,250 160,000 330,000 

R;30;1.00;60 581,726 536,250 278,750 502,500 448,750 

R;30;1.00;80 486,374 403,750 357,500 287,500 587,500 

R;30;1.00;95 341,693 331,250 306,250 373,750 571,250 

R;30;1.00;115 150,554 201,250 210,000 170,000 245,000 
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Table 2. Average values of di. For each instance, the best obtained value is marked 
with f. 

a;f3;j;8 MOLSD PMX OBX UOBX UOBXA2EX 
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Table 3. Average values of di- For each instance, the best obtained value is marked 
with f. 

a;/3;7;<5 MOLSD PMX OBX UOBX UOBXA2EX 
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